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Time-Resolved  Fluorimetry  via  a New  Cross-Correlation  Method 


Abstract : 

A new  instrument  for  the  measurement  of  short  fluorescence  lifetimes 
is  discussed  and  characterized.  The  instrument  is  basically  an  implemen- 
tation of  the  cross-correlation  between  the  excitation  source  and  the  induced 
fluorescence  response.  Such  an  approach,  with  a mode-locked  argon-ion 
laser  as  the  source,  is  shown  to  be  capable  of  measuring  lifetimes  as 
short  as  80ps  with  a precision  of  lOps.  The  instrument  is  simpler  and 
less  costly  than  others  of  similar  capability. 


I.  Introduction 

'WWVAAAAAAA/ \A/W 

Time-resolved  fluorimetry  is  an  extremely  useful  tool  in  a number 
of  branches  of  science.  Not  only  have  f luorimetric  methods  been  instru- 
mental in  understanding  the  photophysics  and  photochemistry  of  atoms  and 
molecules1,  but  they  have  also  proven  useful  for  elucidating  macromolccular 
properties  of  biological  species2’ 3 . Moreover,  time-resolved  fluorimetry 
is  finding  application  in  qualitative  and  quantitative  chemical  analysis4. 

A number  of  instrumental  approaches  which  exhibit  a high  degree  of 
temporal  resolution  have  been  developed  for  the  determination  of  fluores- 
cence decay  curves.  These  include  time-correlated  single-photon  systems5*6, 
phase  shift  methods7*8  and  mode-locked  laser  systems9*10  utilizing  optical 
gating  devices.  Of  these  techniques,  the  most  popular  at  the  present  time 
is  the  time-correlated  single-photon  method.  This  technique  has  very  high 
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sensitivity  but  is  relatively  slow  in  data  acquisition  and  involves  rela- 
tively complex  instrumentation,  making  it  costly  to  implement.  In  contrast, 
phase  shift  methods  are  fast  but  provide  at  a given  source  modulation 
frequency  only  one  measurement  point  on  a decay  curve,  thereby  necessi- 
tating an  assumed  functional  form  for  the  decay  curve.  Finally,  laser- 
based  methods  offer  very  high  temporal  resolution  but  are  costly  and  often 
suffer  from  very  low  signal -to-noisc  ratios9’10. 

In  the  present  paper,  a new  approach  for  the  measurement  of  fluores- 
cence lifetimes  is  introduced  and  characterized.  This  new  approach  involves 
a cross-correlation  of  the  perturbing  light  impulse  with  the  evoked  emission 
response  from  the  fluorescent  species.  In  this  study  a CW  mode- locked 
laser  serves  as  the  excitation  source.  The  cross-correlation  process  is 
implemented  in  an  opt o-clcctvonic  fashion  using  high-speed  photodetectors 
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and  microwave  electronic  components.  The  time-delay  variable  required  by 
the  correlation  process  is  adjusted  by  means  of  an  optical  delay  line. 

This  approach  lias  resulted  in  the  development  of  a time-resolved  f luor inl- 
et er  which  provides  lOps  time  resolution  at  nanomolar  concentration  levels 
while  being  much  less  costly  to  implement  and  simpler  to  operate  than  the 
earlier  approaches  mentioned. 


I 
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II.  Description  of  the  Method 

nAAAAAAAAAA/VVAAAAAAAAnAAAA/VVVV 

The  new  correlation  method  can  best  be  understood  qualitatively  through 
examination  of  the  instrument  designed  to  implement  it.  A schematic  diagram 
of  the  instrument  is  shown  in  Fig.  1. 

' j 

In  the  instrument  a mode-locked  laser  emits  a periodic  train  of  very 
brief  optical  pulses.  A portion  of  this  radiation  is  directed  into  photo-  j 

diodes  PD1  and  PD2  by  beam  splitters  BS1  and  BS2  respectively.  One  of 
these  photodiodes,  PD1,  has  a very  fast  response  time  and  thus  can  faith- 
fully follow  the  temporal  characteristics  of  the  optical  pulse.  The  output  I 

signal  from  PD1  is  sent  to  a multiplying  device,  MULT.  The  portion  of  the 
laser  radiation  transmitted  by  BS1  and  BS2  traverses  an  optical  delay 
line,  consisting  of  corner-cube  reflectors  CC1  and  CC2,  and  impinges 
upon  a fluorescence  sample  in  cell  C.  The  fluorescence  resulting  from 
this  optical  stimulus  is  detected  by  a photomultiplier  tube,  PMT.  The 
fluorescence-induced  photo-current  signal  is  then  sent  to  the  other  input 
of  the  multiplying  device.  The  output  of  this  multiplier  is  a signal 
proportional  to  the  product  of  the  two  input  signals.  The  multiplier 
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output  is  then  sent  to  a signal  processing  device  winch  effectively  per- 
forms a running  time  average  of  this  product  signal. 

Operation  of  the  instrument  shown  in  Pig.  1 can  be  understood  intu- 
itively as  follows.  The  optical  delay  line  formed  by  the  movable  corner - 
cube  reflectors  allows  one  to  change  the  relative  arrival  times  of  the 
optical  pulse  at  the  photodiode,  PD1,  and  the  sample  cell.  By  changing 
this  relative  arrival  time,  a temporal  offset  is  generated  between  the 
photodiode  signal  (which  mimics  the  brief  laser  pulse)  and  the  photomul- 
tiplier output  (which  indicates  the  fluorescence  decay).  Because  the 
electrical  pulse  from  the  photodiode  is  considerably  shorter  in  duration 
than  that  from  the  photomultiplier,  it  serves  as  a sampling  or  gating 
function  at  the  multiplier.  In  essence,  the  output  of  the  multiplier 
is  then  just  a small,  time-sampled  portion  of  the  photomultiplier  signal. 
In  turn,  the  portion  of  that  signal  which  is  sampled  is  determined  by 
the  position  of  the  corner-cube  reflector.  The  output  of  the  multiplier 
then  time  averaged  over  many  such  pulses  at  a number  of  selected  delay 
values  (i.e. , corner  cube  positions)  to  provide  a precise  measurement 
of  the  fluorescence  decay  curve. 

The  process  involved  in  the  foregoing  measurement  is  nothing  more 
than  an  instrumental  computation  of  the  cross-correlation  between  the 
laser  stimulus  and  the  resultant  fluorescence  response.  With  this  fact 
in  mind,  one  can  model  the  behavior  of  the  instrument  mathematically  as 
shown  below. 

In  order  to  determine  a fluorescence  lifetime  from  the  output  of 
the  instrument  of  Fig.  1,  it  proves  useful  first  to  find  the  form  of  the 
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instrument  response  function  itself.  Tlie  instrument  response  function 
is  defined  here  as  that  output  which  the  instrument  would  generate  if 
the  photomultiplier  observed  an  optical  signal  having  the  temporal  char- 
acteristics of  the  laser  itself.  Such  a situation  could  be  created,  for 
example,  by  placing  a scattering  suspension  in  the  sample  cell.  In  this 
case  the  output  of  the  photomultiplier  and  the  photodiode  can  be  described 
as  the  convolution  of  their  respective  impulse  response  functions11  with 
the  function  describing  the  time  dependence  of  the  optical  irradiance. 

The  equations  describing  the  temporal  output  of  the  photodiode  and  the 
photomultiplier,  respectively,  are 


D(t)  = 

fX  hd(y)  S-(t-Y)  dy 

J-oo 

(i) 

r 8 

P(t) 

Aw  vc)  £(t_c)  dc. 

(2) 

where  hd(y)  is  the  impulse  response  of  the  photodiode,  h^(£)  is  the  impulse 
response  of  the  photomultiplier  and  Sl(t)  is  the  time -dependent  laser  irra- 
diance. The  overall  response  for  the  instrument,  I(t),  can  now  be  written 
as  the  cross-correlation  of  the  signals  described  by  Eqs.  (1)  and  (2): 


I(i  ) = A r°°  D (t)  P(t  + T)  dt 

J - oo 


(3) 


In  Eq.  (3),  the  output  of  the  instrument,  I(t),  is  a function  of  the  relative 
delay  parameter,  t.  In  turn,  the  value  of  t is  controlled  by  the  position  of 
the  delaying  elements  of  the  optical  delay  line  mentioned  earlier.  In  the 


above  equation,  A is  a scaling  factor  which  accounts  for  the  gains  of  the 
multiplier  and  the  time-averaging  device.  Combining  Eqs.  (1),  (2),  and 
(3),  one  obtains 

00 

I (x)  = a Iff  hd(Y)  IXpCt)  A(t  - Y)  A Ct  + t - O dCdydt.  (4) 

-00 

Integrating  Eq.  (4)  over  t yields 

00 

I(T)  = A ff  hd(Y)  hp(C)  CU(T  + y - 5)  d?dy  (5) 

-00 

where  is  the  optical  irradiance  autocorrelation  function  of  the  laser 
pulse12. 

Experimentally,  h,(y)  approximates  a Dirac  delta  function.  Using  , 
a I 

this  fact  in  Eq.  (5)  and  integrating  over  y yields  Eq.  (6), 

00 

1 (t)  = A J h (c)  Cu(t  - O dc.  (6) 

-00 

Eq.  (6)  implies  that  the  output  of  the  instrument  when  the  photomultiplier 
observes  scattering  is  simply  the  convolution  of  the  photomultiplier 
impulse  response  with  the  autocorrelation  function  of  the  laser  pulse. 
Significantly,  Eq.  (6)  shows  no  dependence  on  time,  but  only  on  the  optical 
delay  parameter,  x. 

When  the  photomultiplier  observes  fluorescence  rather  than  scattering, 
the  output  of  the  photomultiplier,  P'(t),  is  described  by  the  following 
equation. 

oo 

-00 


E(V)  I’(t  - V)  dv 


(7) 
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In  Eq.  (7),  F(v)  is  the  fluorescence  system  impulse  response  function 
(i.e.,  the  fluorescence  decay  curve)  and  P is  as  described  earlier.  From 
Eqs.  (1),  (2),  (3)  and  (4),  the  instrument  output  for  the  fluorescence 
experiment,  I' CO,  will  be 
00 

I'(T)=A ffff 11(1  CY)  hp(°  F(v)  *(t  ~Y)  l(t  + T " v ■ ° 

-00 

dt  dv  d C dy.  (8) 


As  before  [cf.  Eqs.  (3)  - (6)],  Eq.  (8)  can  be  integrated  over  t and  y, 
assuming  hj(y)  approximates  a delta  function.  This  procedure  yields 


oo 

I'(T)»  hp(0  CU(T  - V - 5)  dvd£.  (9) 

-00 


1 

( 

Comparing  Eq.  (9)  with  Eq.  (6)  (the  latter  being  the  instrument  output 
for  a scattering  experiment  or  for  fluorescence  with  a lifetime  approxi- 
mating zero),  one  observes  that  the  output  for  the  luminescence  experi- 
ment [Eq.  (9)]  is  just  the  convolution  of  the  luminescent  system  impulse 
response  function  [F(v)]  with  the  instrument  response  function  [I(x)], 


i.e. 


00 

I ' (t)  -f  F(V)  I (T  - v)  dv  . (10) 

-00 

Using  the  foregoing  treatment,  one  can  readily  formulate  a method 
for  determining  fluorescence  lifetimes  with  the  new  instrument.  First, 
the  instrument  response,  I(x),  can  be  measured  simply  by  observing  scattered 
laser  radiation..  Next,  I'(x)  can  be  determined  for  the  unknown  luminescing 
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substance.  The  luminescence  lifetime  can  then  he  readily  calculated  by 
deconvolution  of  I'(t)  with  I(t). 

Of  course,  these  deductions  are  limited  to  the  case  of  only  a single 
exciting  pulse  of  laser  radiation.  In  an  actual  experiment  an  "infinite" 
train  of  pulses  is  utilized  to  perform  the  correlation  operation.  Fortu- 
nately, the  previous  treatment  can  be  readily  modified  to  describe  this 
pulse-train  situation. 

The  output  from  a CW  mode- locked  laser  is  a continuous  train  of  optical 
pulses  spaced  apart  where  c is  the  speed  of  light,  n is  the  refractive 
index  of  the  medium  within  the  optical  cavity  and  L is  the  distance  between 
the  two  reflectors  comprising  the  cavity.  Such  a train  of  pulses  can  be 
assumed  to  be  a coherent  pulse  train,  i.e.,  that  all  pulses  have  essentially 
the  same  temporal/amplitude  characteristics.  To  mathematically  express 

l 

this  pulse  train,  the  optical  irradiance  function  described  earlier  must 
be  convolved  with  an  infinite  sequence  of  delta  functions  spaced  equally 
in  time 

/°°  00 

I £(A)  6 C t - nT  - A)  dA  (H) 

n=  -oo 

where  £ is  as  described  previously,  6 is  the  Dirac  delta  function  and  the 
subscript  e denotes  that  the  optical  irradiance  function,  £e(t),  is  now  a 
train  of  the  pulses  described  by  £.  Also  T = i.e.,  T is  the  pulse  spacing 

determined  by  the  laser  cavity. 

To  describe  the  new  instrument  response  function,  Ie(t),  Eq.  (11) 


must  be  inserted  into  F.q.  (4),  which  then  becomes 


where  the  summation  has  been  interchanged  with  integration  and  the  2M 
implies  that  the  time  integration  was  performed  over  time  2MT.  Integrating 
Ilq.  (13)  over  6,  X and  y and  once  again  assuming  that  h^  (y)  -+■  <S(y),  cnc  finds 


I 


00  IKT 


. . _ lim  A 

eTl  = M-x»  2MT 


-oo  -MT 


Z l \ (O  MX)  UP) 

])=-CO  i=»oo  i 

6(t-nT-X-y)  6 (t-iT-P-f.+i) 


The  lim 


d’  ’ydXdPdt. 


arises  because  an  infinite  pulse  train 


is  now  being  correlated13,  and  a time  average  is  being  found  in  Eq.  (12) 
rather  than  a time  integral  as  in  Eq.  (4)'.  Of  course,  to  determine 

Ie(r),  °ne  cannot  afford  to  do  the  infinite  time  average  implied 
by  Eq.  (12).  However,  Ie(x)  will  converge  to  its  infinite  time  average 
value  within  a time  interval  equal  to  a few  cycles  of  the  lowest  frequency 

i 

contained  in  &e(t) . In  the  present  instrument,  this  frequency  is  i 
(the  reciprocal  of  the  mode-locked  pulse  train  period);  accordingly, 
even  the  shortest  practicable  measurement  time  will  produce  a value  of  M ' 
large  enough  to  render  the  approximation  valid.  Conveniently,  t can  assume 
any  value,  because  the  pulse  train  is  essentially  of  infinite  length,  and 
there  is  no  concern  with  end  effects1*. 

Performing  the  time  integration  in  Eq.  (12)  yields 


t , ..  2MA  “ 

‘c(T)  ■ Wf  ,E 

(n-i)=-°° 


hd(y)  h («  UX)  1(8) 


6[ (n-i)T  +X+y-P-C+t]  d^dydpdX 


IC(T)  = ~ z f°°h  CU(T  - C + kT)  dr..  (14) 

k=*°°  J-ao  1 

In  Eq.  (14),  k lias  been  substituted  fox-  (n  - i).  Comparing  Eqs. 

(6)  and  (14),  the  only  difference  is  that  the  latter  has  been  scaled 
by  a factor  of  -jr  and  that  I is  now  a periodic  function. 

Following  the  same  procedure  as  before,  Eq.  (9)  becomes 

oo  00 

I'e(T)  = ~ l JJ F (v)  h (o  CU(T  - v - 5 + kT)  dvdc  (15) 

k=  -oo 

for  the  periodic  pulse  train.  Again  substituting  Eq.  (14)  into  Eq.  (15), 
one  arrives  at  Eq.  (10). 

Several  important  practical  points  derive  from  the  foregoing  treat- 
ment. First,  for  a valid  measurement  there  is  no  need  for  the  change  in 
delay  time  At  to  be  larger  than  T (the  period  of  the  pulse  train)  because 
the  instrument  output  has  this  same  periodicity  in  t.  Second,  the  higher 
the  pulse  repetition  rate, 1/T,  the  greater  the  signal  repetition  rate  and 
the  larger  its  magnitude,  le  or  I'  . Third,  it  is  possible  that  F(t), 
the  fluorescence  decay,  will  be  of  sufficient  duration  that  it  overlaps 
the  next  pulse  of  radiation,  implying  that  information  about  F(t)  will  be 
lost.  In  practical  terms  this  fact  limits  the  length  of  lifetimes  that 
can  be  conveniently  measured. 
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III.  Experimental 

WUWWWWWVWVb 

A.  The  Instrument 

Fig.  1 is  a schematic  diagram  of  the  instrument  used  in  this  work; 
all  components  were  mounted  rigidly  on  a 3’  a 12’  optical  table.  The 
excitation  source  utilized  in  this  instrument  is  a CW  mode-locked  Ar+ 

Laser  (Model  171-06  A r+  laser  and  Model  361  acousto-optic  mode-locker, 
Spectra  Physics,  Mountain  View,  CA) . The  optical  pulse  train  from  this 
laser  has  a repetition  rate  of  approximately  82  MHz  and  a resultant  pulse 
spacing  of  12.2  ns.  Each  pulse  contains  approximately  5nJ  of  optical 
energy  and  has  a full  duration  at  half  maximum  (FDHM)  of  approximately 
150ps  under  the  operating  conditions  utilized  here. 

Immediately  following  the  laser  output  is  a mechanical  chopper  which 
interrupts  the  laser  beam  at  a frequency  of  33Hz  and  a duty  factor  of  12.5 
percent.  The  chopped  beam  intersects  a glass  plate  (BS1)  used  as  a beam 
splitter  to  direct  a small  portion  of  the  optical  energy  into  a fast 
photodiode  (PDl).  The  intensity  of  this  reflected  beam  is  adjusted  with 
a neutral  density  filter  (NDJ).  The  light  is  then  focused  to  a diffraction- 
limited  spot  with  a glass  lens  (LI,  f = 57mm)  and  imaged  onto  PDl . PD1  is 
a specially  constructed  Shottky  photodiode1 1 having  an  impulse  response 
FDHM  of  approximately  50ps.  PDl  is  reverse-biased  with  a type  B (90V) 
battery.  The  alignment  of  this  photodiode  is  quite  critical  because  its 
active  region  is  an  annular  ring  that  is  only  approximately  10pm  wide. 
Illuminating  regions  outside  of  this  active  annulus  causes  charge  separa- 
tion outside  of  the  high-clcctric-f icld  area,  resulting  in  a less-than- 
optimal  time  response.  Accordingly,  PDl  must  he  positioned  exactly  at 
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the  focus  of  lens  LI  and  its  lateral  movement  adjusted  for  host  performance. 
For  positioning  adjustment,  PD1  was  mounted  on  a high-precision  x,  y,  z 
translation  stage  (Model  420-05,  Newport  Research  Corp.,  Fountain  Valley, 
CA). 

Following  BS1  is  a second  glass  plate  BS2.  BS2  deflects  another  small 
fraction  of  the  laser  through  a focusing  glass  lens  (L2,  f = 57nim)  onto 
PD2  (Model  XL56,  Texas  Instruments,  Dallas,  TX) . PD2  is  mounted  directly 
onto  an  SMA  bulkhead  connector  along  with  approximately  . 006pf  of  shunt 
capacitance.  The  reverse-bias  on  PD2  was  supplied  by  two  type  B batteries 
in  series  (180V). 

The  main  laser  beam  is  directed  by  folding  mirrors  Ml  and  M2  toward 
the  variable  optical  delay  line.  Both  are  dielectric-coated  to  reflect 

O 

5145A  radiation  at  45°  angle  of  incidence.  From  M2  the  laser  beam  travels 
through  a variable  neutral  density  filter  (VND)  which  is  used  to  control 
the  optical  power  reaching  the  sample  cell.  The  radiation  then  arrives 
at  the  variable  delay  line,  formed  by  glass  corner-cube  reflectors  CC1 

t 

and  CC2  (Edmund  Scientific,  Barrington,  NJ).  The  return  beam  from  CC1, 
displaced  approximately  3cm  from  the  input  beam,  is  folded  by  the  two 
dielectric-coated  mirrors  M3  and  M4.  Exiting  the  second  corner-cube 
reflector  (CC2) , the  beam  travels  through  glass  lens  L3(f  = 600mm) , neu- 
tral density  filter  ND2  and  into  the  sample  cell.  Corner  cubes  CC1  and 
CC2  were  selected  for  the  optical  delay  elements  over  corner  mirrors 
because  they  greatly  reduce  the  mechanical  tolerances  required  for  the 
delay  line  to  produce  a spatially  stationary  return  (output)  beam. 

Corner  cubes  CC1  and  CC2  arc  mounted  on  optical  holders  which  move 
independently  along  an  optical  rail.  The  rail  is  of  sufficient  length  to 
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allow  approximately  1 meter  of  movement  by  each  corner  cube.  The  resulting 
total  displacement  (Ax  = 2m)  is  sufficient  to  introduce  a relative  time 
delay  of  greater  than  12ns,  which  is  equal  to  the  period  of  both  the  pulse 
train  and  the  cross-correlation  function.  This  optical  rail  is  calibrated 
with  a 1mm  per  division  scale  whicli  allows  one  to  set  easily  any  delay 
time  between  0 - 12ns  with  ±5ps  accuracy. 

The  induced  sample  fluorescence  passes  through  a suitable  interference 
filter  (IF)  and  is  detected  by  an  RCA  C31024  photomultiplier  tube  (PMT) . 
This  PMT  has  a biasing  circuit  which  was  constructed  in  our  laboratory 
and  operates  at  -3500V.  The  biasing  circuit  is  essentially  voltage  divider 
(B)  suggested  by  RCA  for  high-peak-current  operation  of  this  tube.  How- 
ever, the  circuit  was  modified  for  DC  coupling  by  biasing  the  suppressor 
grid  at  approximately  -500V  with  respect  to  ground  and  removing  the 
capacitive  coupling  network. 

The  electrical  signal  from  the  photomultiplier  tube  is  transmitted 
via  RG-58C  coaxial  cable  into  the  RF  port  of  a microwave  mixer  (Model 
ZFM-4,  Mini -Circuits,  New  York,  NY).  This  mixer  is  of  the  double-balanced 
type  with  a bandwidth  of  1.25  GHz  and  an  IF  output  extending  to  DC,  a 

ft 

feature  essential  to  this  application.  Into  the  other  (LO)  port  of  the 
mixer  is  conducted  the  signal  from  the  fast  Shottky  photodiode. 

The  transmission  line  connecting  the  Shottky  photodiode  to  the  mixer 
is  comprised  of  three  main  components.  Connected  to  the  diode  is  a short 
piece  of  RG-58C  coaxial  cable  fitted  with  SMA  cable  connectors.  This 
cable  then  feeds  a constant-impedance,  variable-length  air  dielectric 
coaxial  transmission  line  (Model  874-1, K20I,,  GcnRad,  Concord,  Ma),  whose 
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length  can  be  changed  by  22cm,  corresponding  to  a time  delay  of  0.73ns. 
Following  the  variable-length  line  is  a broadband  (0-8G!lz),  -Cdb  power 
splitter  (Model  874-TPD,  GcnRad,  Concord,  MA) . One  of  the  other  two  ports 
of  this  power  splitter  is  terminated  into  a broadband  50f2  load  while  the 
third  port  goes  to  the  microwave  mixer. 

All  these  components  of  the  transmission  line  between  PD1  and  the 
mixer  are  important  in  the  elimination  of  reflections  of  electrical  pulses 
along  the  line.  Because  of  an  impedance  mismatch  at  the  mixer,  approximately 
10  percent  of  an  incident  pulse  is  reflected  from  the  mixer  port  and  returns 
along  the  transmission  lino  toward  PD1 . In  turn,  the  reverse-biased  photo- 
diode appears  essentially  as  an  open  circuit  and  returns  the  pulse  very 
effectively  toward  the  mixer.  This  process  will  occur  repeatedly  and 
leads  to  a series  of  pulses  rather  than  a single  pulse  arriving  at  the 
mixer.  These  pulses  will  be  separated  in  time  by  an  amount  equal  to  twice 
the  delay  time  of  the  transmission  line;  also,  each  pulse  will  be  approxi- 
mately 10  percent  of  the  amplitude  of  the  previously  reflected  pulse. 

These  pulses  are  largely  eliminated  through  careful  selection  and 
adjustment  of  the  components  in  the  transmission  line.  The  variable  length 
transmission  line  allows  one  to  adjust  the  time  required  for  a pulse  to 
undergo  a complete  up-and-back  reflection  cycle.  When  this  time  is  adjusted 
to  be  exactly  one  pulse-train  period,  each  reflected  pulse  exactly  overlaps 
the  next  original  pulse,  reducing  its  influence  greatly. 

The  power  splitter  inserted  into  the  transmission  line  also  reduces 
the  effect  of  reflected  pulses.  With  the  power  splitter  inserted  into  the 
line,  half  of  the  voltage  of  a pulse  will  be  dissipated  by  the  500  load 
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for  each  traverse  of  the  transmission  line.  Thus  the  amplitude  of  reflected 
pulses  observed  at  the  mixer  port  will  be  lower  in  amplitude  by  4~n  for 
the  n^1  return  pulse. 

The  output  of  the  microwave  mixer  feeds  the  signal  channel  of  a lock- 
in  amplifier  (Model  840  Autoloc,  Kiethlcy  Instruments,  Cleveland,  Oil)  which 
is  terminated  into  a Ml  load.  The  value  of  this  impedance  is  critical, 
because  it  not  only  determines  the  load  presented  to  the  mixer,  but  also 
determines  the  range  of  input  levels  over  which  the  mixer  will  multiply 
properly16 . 

The  reference  signal  for  the  lock-in  amplifier  is  obtained  from  PIN 
photodiode  PD2.  The  output  of  the  lock-in  amplifier  is  then  sent  to  a 
strip-chart  recorder  (Model  SR-204,  Heath  Co.,  Benton  Harbor,  MI). 

Fig.  1 is  consistent  with  the  implementation  of  Eq.  (13)  except  that 

i 

i 

lock-in  detection  rather  than  time  averaging  is  used  on  the  product  signal 
from  the  mixer.  The  equivalence  of  these  two  techniques  can  be  seen  by 
inspecting  Eq.  (12).  Specifically,  time  integration  in  Eq.  (12)  is  of 
the  form 

r 

lim  _1_  / f(t)  dt.  (16) 

T-***  2T  J-x 

Clearly,  the  time  average  or  mean  is  being  found  for  the  time-dependent 
periodic  function  f,  implying  that  the  DC  or  Zero  Hz  component  of  the 
mixer  output  contains  the  information  being  sought.  That  is,  the  magnitude 
of  the  DC  component  of  the  signal  from  the  mixer  output  indicates  the 
magnitude  of  the  cross-correlation  function  at  the  particular  delay  value 
for  which  the  corner  cube  reflectors  are  positioned.  This  DC  signal 


component  can  bo  obtained  as  well  from  a low-pass  filter  as  from  an  inte- 
grator if  measurement  speed  is  not  critical.  To  achieve  an  equivalent 
level  of  performance,  the  RC  time  constant  for  a simple  first -order  low- 
pass  filter  must  be  approximately  8 times  the  integration  time. 

The  lock-in  detection  scheme  is  essentially  a low-pass  filter  approach 
except  that  the  detection  bandpass  has  been  moved  away  from  zero  Hz. 
Chopping  of  the  signal  also  reduces  drift  of  the  photomultiplier  output. 

B.  Operation  of  the  Instrument 

All  optical  and  electronic  delays  were  adjusted  so  the  signals  at 
the  mixer  from  PD1  and  PMT  overlapped  when  one  of  the  corner  cubes  was 
at  its  maximum  delay  and  the  other  slightly  less  than  its  maximum  delay. 
Such  a configuration  allowed  a full  period  of  the  correlation  function  i 
to  be  observed  continuously  across  the  most  prominent  features  of  the 
output.  The  delay,  x = 0,  will  be  defined  here  to  be  when  the  delay  line 
is  set  for  maximum  delay.  In  the  actual  instrument  the  two  respective 
signals  which  overlap  at  x = 0 originate  from  laser  pulses  which  are  3 
periods  apart  or  separated  by  36ns  in  time.  Because  the  correlation 
function  is  periodic,  this  temporal  offset  is  immaterial.  f Tills  case 
corresponds  to  k in  Eqs.  (Id)  and  (15)]. 

The  laser  was  typically  run  at  28  Amps  plasma  current  which  produced 
an  average  radiant  power  output  of  400mW  at  SldSA  under  mode-locked  con- 
ditions. 

Photodiode  PD2,  because  of  its  relatively  slow  response  time,  can  be 
aligned  with  a real-time  500MIIz  oscilloscope  (Model  7A19  vertical  amplifier 


1C) 

with  R7912  transient  digitizer,  Tektronix,  l'nc.,  Beaverton,  OR).  Once 
aligned,  PD2  is  used  to  trigger  the  sampling  oscilloscope  for  observing 
the  output  of  PDl  and  also  supplies  a reference  signal  for  the  lock-in 
amplifier.  PDl  was  aligned  for  optimal  response  with  the  he]])  of  a 75-ps 
rise  time  sampling  oscilloscope  (Model  7T11,  7S11  sampling  plug-ins  with 
S-4  sampling  head  used  in  R7904  mainframe,  Tektronix,  Inc.,  Beaverton,  OR). 

A typical  pulse  from  the  Shottky  photodiode  at  the  mixer,  i.e. , after  the 
power  splitter,  has  a peak  amplitude  of  250mV  and  FDHM  of  200ps  (cf. 

Fig.  2).  With  the  aid  of  a picoammcter  (Model  414S,  Kicthley  Instruments, 
Cleveland,  OH),  the  average  current  from  the  photomultiplier  is  adjusted 
to  approx imately  20yA  with  the  chopper  running.  This  adjustment  is  made  by  chang- 
ing the  variable  neutral  density  wheel  (VND)  and  neutral  density  filter 
(ND2) . 

j 

To  generate  a cross-correlation  function,  the  optical  delay  elements 
were  incrementally  moved  in  50ps  steps  (0.75cm)  from  the  x = 0 to  the 
T = -12.2ns  point.  Thus,  the  strip  chart  recorder  displays  a series  of 
levels,  each  indicating  the  value  of  the  cross-correlation  function,  Fqs. 

(14)  and  (15),  at  a specific  50ps  step.  For  convenience  of  display,  these 
levels  were  re-plotted  to  form  the  smooth  curves  displayed  later  in  this 
paper. 

IV.  Results  and  Discussion 

'WV/WWWt-'WVWWWWWWWV 

A.  Instrument  Response  Measurements 

To  test  the  new  instrument's  time-resolution  capabilities,  laser 


scattering  from  0.22pm  polystyrene  spheres  in  IhO  was  recorded.  Through 
sorb  an  experiment,  the  instrument  response  fcf,  l!q.  (11))  can  he  experi- 
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mentally  determined.  Because  the  optical  irradiance  autocorrelation 
function,  C^,  is  approximately  one-tenth  the  duration  of  the  photomull  i 
plier  impulse  response  lip,  lc  should  be  approximately  equal  to  h . If 
the  assumption  is  made  that  these  functions  arc  Gaussian  in  shape,  the 
FDHM  of  lg  will  be  less  than  one  percent  greater  than  that  of  h itself. 

Fig.  3 compares  the  photomultiplier  response  to  scattering  measured 
(A)  with  the  correlation  approach  to  that  determined  with  a sampling 
oscilloscope  (B)  (x  = 75ps).  The  FDHM  indicated  by  the  cross-correlati on 
instrument  is  1.21  ns  whereas  that  produced  by  the  sampling  oscillo- 
scope trace  is  approximately  1.15ns.  As  expected,  the  cross-correlation 
result  is  very  nearly  the  response  of  the  photomultiplier.  Significantly, 
the  cross-correlation  procedure  faithfully  portrays  details  of  the  photo- 
multiplier response,  including  the  shoulder  on  the  trailing  edge  of  the 
pulse. 

B.  Measurement  of  Fluorescence  Lifetimes 

As  a demonstration  of  the  new  instrument,  several  fluorescence  life- 
times were  determined.  In  these  experiments,  determinations  of  the  instru- 
ment response  and  the  response  to  luminescence  were  alternately  performed. 
In  this  way,  any  variation  in  laser  power  or  drift,  in  electronics  would 
be  apparent.  The  FDHM  of  the  instrument  response  was  found  to  be  repro- 
ducible to  i5ps,  indicating  that  drift  is  not  a significant  problem. 

1.  Convolute-and -Compare  Fitting  Method 

Luminescence  lifetimes  Were  obtained  from  these  data  by  a con- 
v o 1 u t e - a nd - comp  arc  method1.  That  is,  a single  Heaviside  exponential  was 
assumed  to  he  the  impulse  response  of  1 he*  fluorophore  studied.  This 
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function  is  then  convolved  with  the  instrument  response  and  compared  with 
the  measured  fluorescence  response. 

The  computer-based  eonvolute-and-compare  algorithm  was  an  intelligent 
routine  after  that  discussed  by  Zimmerman,  et  al17.  In  this  algoritlvn, 
the  sum  of  the  squares  of  the  residuals  between  the  experimental  and  cal- 
culated curves  is  minimized;  the  amplitude  and  time  constant  of  the  exponen- 
tial impulse  response  are  the  adjustable  parameters.  Briefly,  the  algorithm 
adjusts  these  parameters  by  finding  the  total  differential  of  I ' , the 
fluorescence  response.  The  residuals  for  the  calculated  and  experimental 
curves  at  each  point  are  set  equal  to  this  total  differential  at  each 
point.  The  partial  derivatives  of  the  I'  with  respect  to  the  fit  param- 
eters are  next  calculated.  One  than  has  N equations,  where  N equals  the 
number  of  points  describing  the  curves,  in  two  unknowns  where  the  unknown? 
are  the  differentials  of  the  fit  parameters.  A linear  least-squares  solution 
of  this  overdetermined  system  of  equations  is  performed  to  yield  numbers 
for  the  differentials  of  the  fit  parameters.  These  numbers  are  then  used 
to  correct  the  initial  estimates  of  the  lifetime  and  the  amplitude.  A new 
calculated  curve  is  then  generated  with  the  corrected  parameters  which  in 
turn  generates  new  residuals.  This  procedure  can  be  performed  iteratively 
until  the  sum  of  the  squares  of  the  residuals  reaches  a minimum. 

It  was  found  here  to  be  important  to  adjust  two  additional  variables 
to  obtain  an  acceptable  fit.  first,  a parameter  was  introduced  to  shift 
the  calculated  fluorescence  decay  curve  toward  longer  time  delays.  This 
time  shift  arises  from  the  Stokes  shift  of  the  fluorescence  emission. 

The  higher  energy  of  scattered  laser  photons  yields  photoclcct rons  with 


greater  translational  energies  than  those  produced  by  the  lower  energy 
fluorescence  photons.  The  velocity  vector  of  the  photoolcctrons  collected 
is  predominately  directed  toward  the  first  dynode  in  the  photomultiplier, 
the  net  result  being  that  longer  wavelength  radiation  produces  a longer 
transit  time  for  electrons  to  reach  the  anode18'19.  The  two  wave- 

lengths observed  here  (Stokes  shift  = 85nm)  required  a shift  of  50ps  for 
an  optimal  fit;  in  contrast,  no  time  shift  is  needed  if  AA<30nm. 

Another  parameter  was  introduced  to  enable  the  DC  level  of  the  response 
curves  to  be  adjusted.  This  adjustment  is  necessary  because  the  convolution 
of  a long- lifetime  exponential  with  a train  of  pulses  yields  a curve  whose 
minimum  is  no  longer  that  of  the  pulse  train  but  is  some  greater  DC  level. 
This  effect  would  ordinarily  be  seen  in  the  instrument  output  as  well. 
However,  transformer  couplings  in  the  microwave  mixer  eliminates  any  DC 
level  in  the  present  instrument.  This  DC  adjustment  parameter  was  deter- 
m;ned  to  be  important  only  for  fluorescence  lifetimes  greater  than  a 
nanosecond. 

Optimal  values  for  the  time  shift  and  DC  level  parameters  arc  found 

by  systematically  incrementing  them  through  their  domain  until  a minimum 
* 

of  the  least-squares  sum  is  found. 

2.  Measured  L if climes 

Tig.  4 shows  typical  cross -correlation  plots  generated  i:«  a 
fluorescence  lifetime  determination.  In  each  case  the  instrument  response 
rs  shown  wi th  the  measured  and  calculated  fluorescence  response  curves. 
Immediately  apparent  in  the  instrument  response  plot  of  Fig.  4A  is  a second 
pulse  located  at  vPu.s.  This 


second  hump  is  the  reflected  signal  from 
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the  Shottky  photodiode  which  was  mentioned  earlier.  Fortunately , the 
reflected  pulse  does  not  strongly  affect  the  validity  of  measured  fluores- 
cence lifetimes;  the  lifetime  can  be  retrieved  faithfully  by  the  convolutc- 
and-c ompare  method  as  long  as  all  portions  of  the  detection  system  arc 
responding  linearly.  The  primary  disadvantage  of  the  appearance  of  this 
second  pulse  is  that  it  makes  a log-linear  slope  estimation20  difficult 
even  for  a relatively  long  lifetime. 

Table  I shows  a list  of  the  lifetimes  determined.  A 600nm  interference 
filter,  BW  = 35nm  (Melles  Griot,  Danbury,  CN),  was  used  for  spectral  selection 
at  the  PMT  for  the  first  two  samples  in  Table  I,  whereas  a 546nm  interference 
filter,  BW  = 8.8nm  (Ditric  Optics,  Marlboro,  MA) , was  employed  for  the  last 
two  samples. 

The  two  numbers  for  each  lifetime,  listed  under  LST  SQ  and  DELTA  in 
Table  I,  represent  the  values  calculated  using  two  different  fitting  param- 
eters. The  parameter  LST  SQ  is  the  sum  of  the  squares  of  residuals  or 
normal  least-squares  criterion.  In  contrast,  the  parameter  DELTA  involves 
.ii iiiimi cation  of  the  sum  of  the  absolute  values  of  the  residuals  between  the 
experimental  and  calculated  curves  divided  by  the  integral  of  the  experi- 
mental curve.  Both  of  these  fit  parameters  arc  widely  used25.  Significantly, 
the  results  obtained  using  the  two  fit  parameters  are  not  greatly  different. 

3.  Validity  of  M vjred  Lifetimes 

The  mean  values,  X,  and  a-  -jrd  deviations,  o,  shown  in  Table  I 
are  weighted  values  obtained  from  three  independent  measurements.  These 
weighted  values  were  calculated  using  the  customary  expressions26 
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where  N is  the  number  of  determinations  and  is  the  least-squares  sum 
for  determination  Xj^.  When  the  DELTA  fit  parameter  was  employed,  its 
square  replaced  in  Eqs.  (17)  and  (18). 

Another  informative  parameter  is  the  standard  deviation  of  a single 
lifetime  determination,  0^,  which  is  obtained  using  the  relation26 


a 
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where  M is  again  the  least-squares  sum.  The  derivative  is  found  about 
Jiat  value  of  t for  which  M is  a minimum.  This  quantity  was  calculated 
both  by  deriving  an  explicit  expression  for  aT  and  also  by  numerically 
computing  the  second  difference  of  M about  its  minimum.  For  all  samples, 
was  on  the  order  of  a picosecond,  although  actual  magnitudes  range  from 
slightly  greater  than]  ps  for  Rhodamine  B to  slightly  less  than  1 ps  for  the 
quenched  Fluorescein  sample.  Of  course,  this  number  does  not  define  the 
reliability  of  finding  the  true  lifetime  of  a given  sample  but  instead 
reflects  the  linearity  of  the  detection  system  and  the  amount  of  noise 
in  a given  set  of  data. 


c. 


Instrument  Sensitivity 

The  sensitivity  of  the  present  instrument  and  its  applicability  to 
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dilute  and/or  low-quantum-yield  samples  is  revealed  partly  by  the  operating 
conditions  employed  in  the  preceding  studies.  To  avoid  saturation  of  the 
PMT  by  fluorescence,  the  laser  irradiance  was  ordinarily  attenuated  by 
two  or  three  orders  of  magnitude.  Moreover,  the  lock-in  detection  system 
was  ordinarily  set  to  a relatively  insensitive  scale.  From  this  source 
alone,  three  orders  of  magnitude  greater  sensitivity  could  be  obtained. 

Thus,  lifetimes  of  subnanomolar  concentrations  of  all  samples  examined 
here  should  be  determinable  with  this  instrument.  The  figure  of  merit 
typically  used  to  denote  sensitivity1  is  the  quantum  yield-optical  density 
product  required  to  generate  a detection-limit  signal  for  a sample  of  1 cm 
path  length.  Under  the  conditions  used  in  this  investigation,  this  figure 
of  merit  is  calculated  to  be  between  10”4  and  10~s.  Time-correlated  single- 
photon instruments  (the  most  sensitive  technique)  have  figures  of  merit 
near  10  7 1 . 

D.  Practicality  of  the  New  Method 

' The  new  instrument  is  both  simpler  and  less  expensive  than  others  of 
similar  capability.  Lifetimes  as  short  as  80ps  have  been  shown  to  be 
measurable;  lifetimes  as  long  as  50ns  should  also  be  measurable,  this 
limitation  being  due  to  the  periodicity  of  the  pulse  train.  Lifetime 
measurements  with  this  instrument  also  exhibit  high  precision,  being  as 
good  as  ±10ps  for  the  shorter  lifetimes. 

The  new  instrument  lends  itself  well  to  automation.  A microprocessor 
could  control  the  scanning  of  the  optical  time-delay  elements  while  d ig  i t i 7.ing 
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the  output  of  the  lock-in  amplifier.  The  instrument  also  avails  itself 
to  signal  averaging  if  automation  is  introduced.  Calculating  the  maximum 
rate  at  which  the  delay  elements  can  be  linearly  scanned  such  that,  the 
bandwidth  of  the  changing  "DC"  component  from  the  mixer  falls  within  the 
-3db  bandwidth  of  the  lock-in  amplifier  low-pass  filter  allows  one  to 
determine  the  time  needed  to  do  a complete  t scan.  With  the  300ms  time 
constant  typically  used  for  the  lock-in  amplifier,  one  finds  that  a complete 
scan  from  0 to  12.2ns  can  be  performed  in  less  than  6 sec.  Thus  signal 
averaging  could  be  employed  for  improving  signal-to-noise  ratios. 
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Tabic  I.  Fluorescence  Lifetime  Results 
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